Induced Random Fields in the LiHo x Yi_ 2; F4 Quantum Ising Magnet in a Transverse 
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The LiHo^Yi_a:F4 magnetic material in a transverse magnetic field B x x perpendicular to the 
Ising spin direction has long been used to study tunable quantum phase transitions in a random 
disordered system. We show that the B x — induced magnetization along the x direction, combined 
with the local random dilution-induced destruction of crystalline symmetries, generates, via the 
predominant dipolar interactions between Ho 3+ ions, random fields along the Ising z direction. This 
identifies LiHo I Yi_ a; F4 in B x as a new random field Ising system. The random fields explain the 
rapid decrease of the critical temperature in the diluted ferromagnetic regime and the smearing of the 
nonlinear susceptibility at the spin glass transition with increasing B x , and render the B x — induced 
quantum criticality in LiHo^Yi_ a: F4 likely inaccessible. 
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Quantum phase transitions (QPTs) occur near ab- 
solute zero temperature and are driven by quantum 
mechanical fluctuations associated with the Heisenberg 
uncertainty principle and not by thermal fluctuations 
as in classical phase transitions ll|, The trans- 

verse field Ising model (TFIM) [J \M with Hamilto- 



nian Htfim 



"5 J ij a l°] ~ r *f» where a i 



(ji = x,y,z) are Pauli matrices, is the simplest theo- 
retical model that exhibits a QPT [H, 0, H . The field T 
transverse to the Ising z direction causes quantum tun- 
neling between the spin-up and spin-down eigenstates of 
er? . These spin fluctuations decrease the critical temper- 
ature T c at which the spins develop either conventional 
long-range order or, for random ferromagnetic and anti- 
ferromgnetic Jij, a spin glass (SG) state with randomly 
frozen spins below a, T g . At a critical field r c , T c or T g 
vanishes, and a quantum phase transition between either 
a long-range ordered or SG state and a quantum param- 
agnet (PM) insues. 

The phenomenology of both the disorder-free and ran- 
dom TFIM has been extensively investigated in the 
LiHo 2; Yi_ 2 ;F4 magnet with a magnetic field B x applied 
transverse to the Ho 3+ Ising spin direction which 
is parallel to the c-axis of the body-centered tetragonal 
structure of LiHo x Yi_ :E F4 Crystal field effects give 
an Ising ground state doublet, |$q ), and a first excited 
singlet, |<I> C ), at approximately 9 K above the ground 
doublet [9j. For x = 1, LIH0F4 is a dipolar-coupled fer- 
romagnet (FM) with T c = 1.53 K Random dis- 

order is generated by replacing the magnetic Ho 3+ ions 
by non-magnetic Y 3+ . Quantum mechanical (spin flip) 
fluctuations are induced by B x which admixes l^o) with 
|$^), splitting the latter, hence producing an effective 
TFIM with T = T(B X ) (r oc B x for small B x ) [n|. 

Two experimental puzzles pertaining to the effect of 
B x on the FM (0.25 < x < 1.0) and the SG (x < 0.25) 



phases of LiHo x Yi_ 2 ,F4 have long been known. Firstly, 
in the FM regime, while the mean-field argument that 
T c (x) oc x for the PM to FM transition is well satisfied 
for 0.25 < x < 1.0 8], the rate at which T C (B X ) is de- 
pressed by B x becomes progressively faster than mean- 
field theory predicts as x is reduced [UJ]. This implies 
that, compared with the energy scale for FM order set by 
T C {B X = 0), B x becomes ever more efficient at destroy- 
ing FM order the lower x is ll|. Secondly, for B x = 0, 
LiHoo.i67Yo.833F4 displays a conventional SG transition, 
with a nonlinear magnetic susceptibility X3 diverging at 



T g as xaiT) cx {T-T g )-"i [12j. However, X a(T) becomes 
less singular as B x is increased from B x = 0, with no in- 
dication that a QPT between PM and SG states occurs 
as T — > 0, Uil]. It has recently been suggested that 
for dipole-coupled Ho 3+ ions nonzero B x generates both 
longitudinal (along the Ising z direction) T4] and trans- 
verse [15l |16| random fields (RFs) that either renormalize 
the critical transverse field [I5j], or even destroy the 
SG transition (l6j . In this paper we examine the quanti- 
tative merit of this hypothesis by comparing results from 
numerical and analytical calculations with experimental 
results on LiHo 2; Yi_ 2: F4. We obtain compelling evidence 
that RFs are manifestly at play in LiHo 2; Yi_ ;E F4 and ex- 
plain the above two long-standing paradoxes. 

We first show that the low-energy effective theory of 
LiHo^Yi-^ for x < 1 and B x > is a TFIM with 
additional B x — induced RFs. We start with Hamilto- 
nian Ti. = Ho + 7idip expressed in terms of the angular 
momentum operator J of Ho 3+ (J = 8, L = 6, S = 2). 
The single ion part, H Q = Z)jWcf(Jj) + Hz(Ji)], 
consists of the crystal field Hamiltonian, HrfjJj), of 
Ho 3+ in the LiHo 2; Yi_ a ;F4 environment 0, [juj, and 
the Zeeman field term, Hz = —g^{Ji • B), with 
B the magnetic field, g = 5/4 is the Ho 3+ Lande 
factor and ^ is the Bohr magneton. The inter- 
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actions between ions are dominated by long-range 
magnetic dipolar interactions 7idip @, Gj}' Ti-dip = 
(<? 2 M|/2) E(i,j) [ J * • J j - 3(J< • r^Jj ■ rij)rr. 2 ] r„~ 3 , 
where are the crystalline positions occupied either by 
a magnetic Ho 3+ ion (e, = 1) or a non-magnetic Y 3+ ion 
(ej = 0), and ry = |rj — r»| is the inter-ion distance. 

The two lowest energy eigenstates \^f(B x )) and 
\^~(B X )) of 7io are sufficiently below |$ e ) that the 
latter can be ignored at temperatures near and be- 
low T c (x = 1) - 1.5 K [13]. This allows us to re- 
cast TL in terms of an effective Hamiltonian, 7i c ff, 
with S = 1/2 pseudo spin operators that act in 
the restricted low-energy subspace spanned by the 
Hi l^rt 8 ^)) fa = ±) eigenstates of TL . In this 
subspace, T(B X ) = (l/2)[(*+|H |* + > - (*-|Wo|*~>]- 
The projected Jj* (/i = x, y, z) operator is written as: 
J? = Y,u C A B x)°i + C ^t 0. The |+) and |-) 
eigenstates of of are written in terms of \^f) such that 
J[- = C zz af [k|. For B x = 0, only C zz ^ and decreases 
slightly with increasing B x , while the other C I1V (B X ) pa- 
rameters and T(B X ) increase with B x , starting from zero 
at B x = 0. By straightforward manipulations replacing 
Jf in terms of Ei/ ^W^f in TLdip, one finds that the 
terms with largest C )1V {B X ) in TL e e, including the trans- 



verse field term Taf , are: 



n, 



elf 



c. 
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where L 



ah 



[Sab 
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3 . We see from Eq. [1] 



that a longitudinal random field (RF) term oc of along 
z and an off-diagonal (bilinear) oc of coupling are in- 
duced by B x ^ (C x q — C xx = for B x — 0). For pure 
LiHoF4 (a; = 1, all tj = 1), lattice symmetries enforces 
J2j = 0, causing the term linear in of and the Boltz- 
mann thermal average (of cr|) to vanish for i?^ > 0. How- 
ever, for x < 1 and > 0, a RF oc af emerges. With the 
time-reversal symmetry broken by B x ((Jf) > 0), the bi- 
linear of of also provides a "mean-field" contribution to 
the longitudinal RFs, (of )aj, as well as transverse RFs, 
(of >of . We find that C xx /C xQ < 1 for all B x > so that 
the leading correction to -Htfim is indeed a (correlated) 



RF term, ~ E, ^of , with oc C x0 C zz ^ ; . , 

We now investigate the effect of the RFs on the PM 
to FM transition. To do so, we make a mean-field 
(MF) approximation to TL and consider the one-particle 
MF Hamiltonian, TL^ , for an arbitrary site r$ occu- 
pied by a Ho 3 + moment: Hf v = H ci - gm{hf Y ■ 
Jj) — g^JfB x . h.^ F is the MF acting on magnetic mo- 



ment Ji, h 



MF 



1 y ' 

with the self-consistent MF equation Mi — (Jf) with 
M, - Tr(pM F J J )/Tr(p M F) and ^mf = exp(-/3iJf F ) 
with /3 = \/(k-B,T). The crystal field Hamiltonian H c { for 
Ho 3+ expressed in terms of the components Jj* is taken 
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FIG. 1: T C (B X ) from a finite lattice mean-field calculation 
(see text). For x = 0.7 and x — 0.5, the filled symbols show 
T C (B X ) when only the Lff dipolar terms that couple z com- 
poments of Ji are kept. The open symbols show the increased 
rate of depression of T C (B X ) when the off-diagonal dipolar Lf/ 
couplings are included. For the pure x = 1 case, T C (B X ) is the 
same for both models since the lattice symmetries eliminate 
the internal local random fields along z. 



from the 4 crystal field parameter (4CFP) H c { of Ref. [9(. 
Slightly different choices of H. c f [a, d° not qualita- 
tively affect the results. We diagonalize 7if IF for each 
i and update the MiS using the newly obtained eigen- 
states. We continue iterating until convergence is reached 
at the n'th iteration, defined by the convergence criterion 

Ei( M f n) - M ( l n - 1) ) 2 /E t (Mt ) ) 2 < 1(T 6 . To simplify 
the calculations and speed up the convergence we only 
keep the diagonal Lf? and the off-diagonal Lf? terms in 
Hdip since no qualitatively new physics arises when keep- 
ing the other diagonal (Lf? and L\J) and off-diagonal 



(if 



, ( . „ y ) terms in Hdi P - For a given B^, we com- 
pute the temperature dependence of the magnetization, 
M z = [Ei M ih, along the I direction Here [. . .] d 

signifies an average over the bimodal lattice occupancy 
probability distribution, P(ei) = xS(ei — 1) + (1 — x)S(ei). 
The transition temperature T C (B X ) is determined by the 
temperature at which M Z (T) sharply rises [lil ]. We con- 
sider a system of linear size L — 4, with total number of 
sites iVo = 4L 3 with N = xTVq sites occupied by Ho 3+ , 
and perform disorder averages over 50 different diluted 
samples. We use the Ewald summation method to define 
infinite-range dipolar interactions. For B x = we find 
that, as found experimentally, the decrease of T c while 
reducing x is proportional to x d, [ll|. To investigate 
the role of RFs, we compare T C (B X ) obtained when both 
the Lf? and Iff terms in 7idi P are kept (open symbols) 
with the T C (B X ) found when only the Lff term is re- 
tained (filled symbols). Figure 1 shows that T C (B X ) is 
depressed for small B x faster when the off-diagonal Iff 
term in 7idip is present than when only the I sing iff 
term is kept. Also, as found experimentally ll|, the 
rate at which T C (B X ) is depressed by B x increases as x 
is lowered. This provides strong evidence that the ex- 
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perimental observation is due to B x — induced RFs whose 
variance increases as x decreases or B x increases. 

We now consider the role of RFs at the SG tran- 
sition. While the nonlinear susceptibility X3 diverges 
at T g ~ 0.13 K in L1H00.167Y0.833F4 when B x = 0, 
Xz(B x ,T) becomes progressively smeared as B x is turned 
on (see top inset of Fig. 3) @, [l3f . It has been sug- 
gested that random off-diagonal dipolar couplings de- 
stroy the SG transition when B x > [16|. Here, we 
take a more pragmatic approach and ask whether the be- 
havior of Xs{Bx > 0, T) in SG samples of LiHo :r Yi_ :r F4 
can indeed be interpreted in terms of induced RFs. To 
investigate this question on X3(B Xl T), we introduce a 
mean- field variant of Ti c s in Eq. (1) that preserves the 
crucial physics therein. Our model is a generalization 
of the Sherrington-Kirkpatrick transverse field Ising SG 
model [B|, [3, [20, but with additional off-diagonal and 
RF interactions similar to that in TL c q: 



n 



\ J ij a i a j 
(i,j) 



(2) 



For simplicity, we take infinite-ranged Jy and Kij given 
by independent Gaussian distributions of zero mean and 
variance J 2 and K 2 , respectively, /if is a Gaussian RF 
with zero mean and variance A 2 . Model (2), but with 
— 0, was previously used to calculate X3 m quadrupo- 



lar glasses [19( , which also possess internal RFs [22| . 

We follow the procedure of Ref. [21|, employing the 
imaginary time formalism and the replica trick to de- 
rive the (replicated) free energy of the system. To fur- 
ther simplify the calculations, we make a static approx- 
imation for the replica-symmetric solution in the PM 
phase. This allows us to derive self-consistent equa- 
tions for the a components of the magnetization, M a , 
and spin-glass order parameters Q a (a = x, z). We find 
M a = (1/2tt) f x oo dxdze- ( - x2+z2 '>/ 2 [(H a /H)tanh(j3H)]P 
with p = 1, H z = (h z + z^/J 2 Q z + (K 2 /2)Q x + A 2 ), 
H x = (r + XyJ(K 2 /2)Q z ) and H 2 = H 2 + H 2 . The 
self-consistent equations for Q x and Q z are obtained by 
replacing M a — > Q a above and setting p = 2. X3 i s 
obtained from X3 = \d 3 M z / {dhl) 3 , and by solving nu- 
merically the resulting four coupled self-consistent equa- 
tions. Figures (2a) and (2b) show the T dependence of 
X3 for various temperatures, T, in models either with 
only hf random fields (Fig. (2a)) or with only random 
off-diagonal K%j couplings (Fig. (2b)). 

Comparing Fig. (2) with the experimental X^i^x) 
(top inset of Fig. (3)), one finds that the dependence 
of X3 upon T in Fig. (2a) does not show a decrease in 
magnitude as T is decreased. Also, while X3 shows a 
decreasing amplitude with decreasing T in Fig. (2b), it 
does not reveal a rapid sharpening as T is increased. It 
turns out that the key physics ingredient missing in these 
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FIG. 2: (a), top panel: \3 vs r with K = 0. (b), bot- 
tom panel: X3 vs F with A = 0. The change of linestyle at 
low r indicates the limit of stability of the replica symmetric 
paramagnetic solution, as determined following the standard 
procedure [2H ]. 



calculations is the underlying microscopic dependence of 
J, A and K upon B x via the C^ V {B X ) transformation 
coefficients. Physically, the built-in B x dependence 
of the Cfxv insures that the B x scale at which X3 is 
quenched by RFs is not trivially tied to the scale at 
which the T = SG-PM crossover occurs, as it is in H 
with Jij, K^, T and h\ independent of B x (c.f Fig. (2)). 
The widths J(B X ), A(B X ) and K(B X ) are obtained by 
calculating the disorder average of the first, third and 
fourth lattice sums in Eq. (1), respectively. We have K 2 
= A(g^nC zz {B x )C xx {B x )] 2 [^ {l . 3] ^(^ z ) 2 ]d and 
A 2 = (g^nC zz {B x )C M {B x )] 2 [^- Q ei(E jyi e,-Lf/) 2 ] d . 
To make further contact between calculations and ex- 
perimental data, we note that hyperfine interactions, 
which are important in Ho-based materials 0, fltj| , lead 
to a renormalization of the critical B x , B x , when T c or 
T g is less than the hyperfine energy scale [lj| EH- To 
obtain a relation between T and B x in li., which does 
not incorporate hyperfine effects, we set T{B X ) / J(B X ) = 
1.05(B X /B^)°- 35 , where B% = 1.2 T is the experimental 
zero temperature critical field This ansatz for 

Y(B X ) is obtained by matching the critical temperature 
T g (B x ) of LiHozYi.^ with the T g (T) of H. For 
the former we use T g (B x ) = TJ0)[1 - {B x /B c x f} 
(4> ~ 1.7) as found experimentally [13j. For Ti., we find 
T g( r ) ~ T s(0)[l - a (r/T g (0)) V '] (a « 0.79 and ip « 4.82) 
by fitting T g (T) vs T with K = A = 0. To incorporate 
the role of hyperfine effects on hf and K^, we rescale 
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T 25 (mK) 
O 66 (mK) 
A 82 (mK) 
O 98 (mK) 
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FIG. 3: X3 hc ° vs B x for different temperatures given by the 
model Eq. (2). The top inset shows X3 XP from Ref - Hi- The 
parameters K and A are computed using x — 0.167 and 
rescaled by n — 0.15 are shown in the bottom inset (see text). 



K and A calculated from their microscopic origin in Eq. 
(1) by a scale factor r\ = 0.15 that positions the peak of 
xf c °{B x ) at B x ~ 5 kOe for T g (B x )/T g (0) = 0.75. With 
an estimate of r\ available, one could then determine a 
parametric T(B X ) such that, for a given T/T g (B x = 0), 
the experimental Xg Xp and theoretical X3 he ° peak at the 
same B x value for all T g (B x )/T g (0). Such procedure 
gives results qualitatively very similar to those in Fig. 3. 

-^thoo reproduces the overall trend of Xg Xp This is ev- 
idence that RFs are indeed at play when B x > and 
at the origin of the experimental Y3 (B x , T) behaviors in 
SG samples of LiHoj;Yi_ x F4 A noticeable differ- 

ence between x 3 hco and is that, for fixed T/T g (0), 

X3 1100 collapses more rapidly and at smaller B x than X3 Xp ■ 
This is likely further manifestation of the renormaliza- 
tion effects of the RFs and random off-diagonal couplings 
caused by the aforementioned hyperfine interactions. We 
will report elsewhere results addressing this hyperfine 
renormalization of Tt c g and its role on X3(B X ,T). 

In conclusion, by comparing numerical and analytical 
results with experimental data, we have obtained com- 
pelling evidence that induced random fields (RFs) are 
indeed at play and "observed" in LiHo x Yi_ x F4. As a re- 
sult, LiHo x Yi_ x F4 in a transverse field (TF) is identified 
as a new RF Ising system. As found in other RF sys- 
tems, we expect that the nontrivial fixed point of the 
theory (7i c ff) is controlled by a fluctuationless classi- 
cal zero temperature fixed point. In particular, this is 
what occurs in a TFIM plus RFs (H with K %i = 0) [13]. 
Hence, it would therefore seem that quantum criticality is 
most likely innaccessible in ferromagnetic LiHo a; Yi_ :E F4 
samples. For the Ising spin glass (SG) model with RFs 
along z, recent studies suggest that there is no Almeida- 
Thouless line 

El 

and no thermodynamic SG transi- 



tion 24|. Hence, from the arguments above leading 
to TC e B, B x — induced SG to PM quantum criticality in 
LiHo :E Yi_ 2 ;F4 would also appear likely inexistent. The 
presence of B x — induced RFs and the quenching of quan- 
tum criticality is presumably the reason why, unlike in 
quantum Monte Carlo simulations of TF Ising SG mod- 
els [2^], Griffiths-McCoy singularities 25 1 have not been 
reported in LiHo^Yi-a^ [aS]- While quantum critical- 
ity seems unlikely for any B x > and x < 1, a new set 
of interesting questions has nevertheless arisen: do FM 
samples of LiHo x Yi_ x F4 with B x > indeed display 
classical RF criticality and all the fascinating phenom- 
ena of the RF Ising model [26]? Given the scarcity of 
real RF Ising materials [13] , the identification of another 
such system opens new avenues for future theoretical and 
experimental investigations. 
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